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ABSTRACT 


A theoretical and experimental investigation of the stress dis- 
tribution and deformation of twisted blades under combined tension and 
torsion was conducted by LCDR Lawrence M. Hogan, USN, at the Aeronau- 
tical Engineering Department of the U. S. Naval Postgraduate School, 
Monterey, California. Theoretical expressions were developed to relate 
the deformation of a twisted "blade" of rectangular cross section and 
constant rate of twist to an applied load and moment. These expressions 
were found to be linear for the rates of twist investigated and depended 
only on the blade material and geometry. The theoretical strain distri- 
bution was found to be essentially parabolic across the chord, with the 
shape determined by the initial twist rate and the angular deformation, 

The experimental investigation was conducted with three specimens 
of constant cross sectional aspect ratio but different rates of twist. 
Correlation with theory was found to be excellent for the chordwise dis- 
tribution of strain and generally good for the linear relationships 
linking the deformation with the applied loads, Recommendations were 


made for areas of further investigation. 
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I» Introduction. 


An analysis of the stress distribution and the distortion of twisted 
blades subjected to combined tensile and torsional loads has several fun- 
damental engineering applications. Propeller and turbine blades are two 
examples of the special case in which the blade is under tension but es- 
sentially free of torsion. High rotational velocities place these blades 
in tension due to centrifugal force while only an insignificant amount of 
torsion results from the aerodynamic loads on the blade. The need to 
predict the resultant stress distribution and accompanying distortion 
of these blades is evident. 

This paper is an attempt to evaluate the applicability of a simple 
theoretical approach to this problem. Specifically, the analysis involves 
a twisted "blade" of constant rectangular cross section and constant rate 
of twist, subjected to known applied torsional and tensile loads. The 
predicted results for three different rates of twist are compared with 
experimental results obtained from three test specimens. 

The strain at any chordwise station is determined from an analysis 
of the geometry of the twisted blade. All deformations are restricted 
to the elastic region. The central chordline of the blade passes per- 
pendicularly through the longitudinal axis. A key assumption is that 
this chordline, originally straight, remains straight in the deformed 
state. Thus the displacement of a chordline a small distance from a 
reference chordline is expressed as a linear displacement in the longi- 
tudinal direction and a rotation about the blade axis. The strain at 


any chordwise station is established as being dependent on ititial rate 


of twist, on the longitudinal strain at the centerline, and on the 
angular deformation. 

The basis for the theoretical analvsis can best be described as 
modified membrane theory. Expressions are obtained for an applied ten- 
sile force and for an applied moment in terms of membrane stress and 
the associated geometry of the blade. The stresses are then expressed 
in terms of the previously determined strain relationships by means of 
Hooke's Law. This allows the applied force and moment to be expressed 
in a form linearly dependent on the longitudinal strain at the center- 
line, e and on the non-dimensional angular deformation, 5 c These 
solutions are modified by the torsional relationships for a tvisted pris- 
matical bar of rectangular cross section. The resultant linear expres- 
sions are then inverted to solve for e, and X in terms of the applied 
tensile and torsional loads, 

Test data was obtained to compare the actual relationships with 
those predicted bv theorv for two special cases of applied load and 
moment: l) An applied tensile load with zero applied moment and 2) 


an applied moment with zero applied tensile load. 
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23 Analysis. 


A diagram of the model to be analyzed is shown in Fig. 1. The 
model is considered to be an axisymmetric helical blade machined from 
a circular rod. The angle of twist, represented by 9 in view (a), varies 
uniformly over the length, L, giving ə constant, The edge of the 
blade describes a helix, the projection of which is a sine wave as is 
evident in view (b). Any cross section of the blade normal to the 
z-axis is essentially a rectangle of length 2a and width..h. 

The helix angle, f, describes the angular difference between the 
tangent to the helix at any radius, r, and a line parallel to the z-axis, 
as shown in view (c). Thus Q varies linearly with r from zero at r = 0 
to سے‎ at r = a. 

As a first approximation in the analysis, the blade thickness is 
assumed very small so that it approaches a pure membrane. The blade is 
therefore considered to have zero bending stiffness and the applied load 
must be taken by tensile stresses in the center surface of the blade. 

An inspection of the applied load, P and T, shows the stress in the 
radial direction to be zero. Therefore, the significant stress for this 
analysis is that represented by gu the stress directed along the tan- 
gent to the helix. Considering then the element of cross section surface 
on which S acts, tdr, the elemental membrane tensile force is repre- 
sented by T tdr, view (c). The component of the tensile force in the 

z direction, g „tdr cos Ø, integrated over the chord must be equal to 
any applied tensile load, P. The component of the tensile force perpen- 


dicular to the z-axis, T tdr sin Ø, introduces a moment about the 
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centerline since it acts at a distance r from the axis. Thus the inte- 
grated moment over the chord must be equated to the applied moment, T. 
These relationships are expressed in eqns. (1) and (2). 

< 


g, tdr cos 67 (1) 
— 


a 
Јоха sin @ (2) 
-Q 


The assumption of a rectangular cross section is not exact in actu- 


FJ 
u" 


= 
it 


alitv due to the machining process. The blades herein analvzed were 
machined from round stock, The stock was mounted in a manner so as to 
pass spirally by a cutting bit, as illustrated in Fig. 2. If this 
cutting bit were of infinitesimally small radius, it would act only at 
a. particular z-station and an exact rectangle would be formed in the 
cross section normal to the z-axis (neglecting the arc on the ends). 

The finite radius of this bit caused it to machine on adjacent z-stations 
as well, introducing a variation in the thickness across the cross sec- 


tion. Thus the thickness t can be found from Fig. 2 to be 
t = h cos 6 - 2R(1 - cos 6) (3) 


where R is the radius of the cutting bit, 
Rearranging and defining the non-dimensional cutter bit radius, 


7” , as R/h gives 
t = h [cos 0 - 2 ۶٢ )1 - cos 2 (4) 
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Fig. 3 shows an elemental length of blade, dz. Line OP is the 
helical path at radius r over the elemental length. The circumferential 
displacement over the length dz is seen to be rdO. Thus the helix angle, 


Ø, is defined by 


0 = asa rdə 


= - (5) 


The non-dimensional rate of twist, K, and the non-dimensional radius, 
f are defined as 
de 


K= a q (6) 


۳ 7 (7) 
Thus eqn. (5) becomes 

f'e tan ` Kf (8) 
giving 

2. тк (9) 

aun Ge nur (10) 

na : )11( 
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2 / + K2 2 
l =- cos ó < / - =—— =- 一 ا‎ (12) 


The radical in the numerator of eqn. (12) can be expanded by the 


binomial theorem to give 


4 
Visawtg? = / 2 22 - گهای‎ + >> )13( 
Dropping higher order terms there is 


4 kif | 
1 - cos 0 = Vig (14) 


Thus eqn. (4) becomes 


I - + kK3 £ 


م - 


(15) 


Substituting these expressions back into eqns, (1) and (2) and 


changing the variable of integration from dr to d£ , produces 


j _ 2 م‎ 2 
۳ 2ah İq, 2070 (16) 
T. nl a EC IEF) ع0"‎ ۶ 2 ٦ dg (17) 


Lf ۶ 22 ”” 7 


Nov an expression for ə. in terms of deformations vill be de- 
rived. Referring once again to Fig, 3, an important simplifying assump- 
tion is made which claims that the radius SP, originally a straight line 
normal to the z-axis, remains a straight line normal to the z-axis in 


the deformed state, S'P', 
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The unstrained length of an element along the tangent is 


c ae ٠ ala V /+ 2 2 (19) 


cos ğ 


The strained condition is defined by a change in length parallel 
to the z-axis, du, and an angular change, dY . The length along the 


tangent to the helix in the strained condition, OP', is found as follows: 


V (dz ita)” +[r(d@+ dy) |* 
oP' = dz V(/+u’)* +[r(o’4 (> 


dz Ver un? 7 — « *2,)2 


ОР” 


۱ 


GEN 


4 
Replacing u' and Vol by the non-dimensional values e, and X/K, respec- 


tivelv, and expanding gives 


ОР" = dz V/ z2e,+ e, + K2g2( / + 2 2 + Х/,г) 
ов" = dz V (/en€2) +[e,(2+e,) + Zí g2(2 6 7«)| 


a) + Xk$”(272 
de ETA 


0۳۱ )20( 


ti 





The strain along the tangent is now expressed bv 6٠ 


Mo 2 
— OP 
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e (2+e,) + Z w £ (2 + ук) 
. 7277 “ 


The radical may be expressed as the binomial expansion of (1 4 p)? 


and substituted as follows: 


(0 
ct 
4 
— 
— 
+ 
Ni 
T 
l 
беј-- 
n 
to 
+ 
— 
| 


e = ی‎ FEF - + F a ہے‎ (22) 


Eqn. (22) is therefore the strain along the tangent to the helix 
at anv radius, f . An order of magnitude analysis of the F term is 
now desired for simplifying purposes. e, will be very small compared 
to unity, as will be X/K if "large" initial twist rates are consid- 
ered; i.e., if the ratio of angular deformation to initial twist is 


small. Neglecting these terms, F in eqn. (22) becomes 


e, + UK E* 


jap = ez‏ اس 


The denominator of Fi is of the order of magnitude unity while the 
numerator is of much smaller order. It is therefore possible to neglect 
powers of Fi and the expression for e, becomes 


2 
Бәј ии (24) 


orun ER 


Eqn, (24) is the simplified expression for the strain along the tangent 


at radius E . It can be rearranged to produce a linear relationship 
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between e. (1 + vi 7 ə and £ ° of the form, y = b + mx as follows: 
2 2 2 
e, (L-KE J) =e, + (K XX) f (25) 


Since the strains are limited to the elastic region, Hooke's Law 


may be applied to give 


Су = Ев, 


Co + Kx g* 
Z+ = — 26 
¿= E 7 (26) 


Eqn. (26) substituted into eqns. (16) and (17) produces 


/ 
(êg SLE e =) 
| ردير‎ A 


Tuğ 2 ع با‎ EN ۲222۸-7۱2 2۶21 Z 28 
= Q ~ E — f ( ) 


)27( ع اه 2 P=‏ 


Non-dimensionalizing these expressions there is 


p 








8 سس 
ZahbE‏ 
eK? l-mkijifl‏ . 
P* - (29)‏ 
ләри“ SEL f‏ | 
| ات TAE‏ 
E‏ چا گم 2 
f‏ 
i=. e f m$ - или, re? əə (30)‏ 
msc ۶2 E E Ug ۳‏ 
f‏ ” 
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Eqns, (29) and (30) are the fundamental theoretical equations to 
be solved. The various integrals are constants determined by the non- 


dimensional rate of twist, K, and are redefined as MK coefficients 


so that 
рх = < i Co + ex 2 E (31) 
T* = << 2 1 Es t eX 22 x (32) 


where X19 = Cé A1: 


These expressions can be inverted to solve for e, and X. 


کی 
e = ————OnƏQ 9?” +‏ 


-— 2 2 
= T” (33) 
0 Ku Aza 7 e 


9, وم مه‎ Teu 


^ 


— ot OC | t 
A фу ии ہیں ے‫‎ 
ex, 22 = əd iz e eX 22. 7 e 2. 


The coefficients of P* and T* in these equations are also constants and 


are redefined as B coefficients. Hence 


иә En p' 7 m e" (35) 


H 


x DA p' T Boo qus (36) 


where Bu: B. 


Solving for the e coefficients in terms of 8 coefficients gives 
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Bea 


cx E E ua s 
11 
852: SN 6.2 Ba 
Ce = جرم‎ 
ҹ Bu Boa 5 Bia Bas 
- Ba, 
en 


gu Вам > Bra Ba 
Bu 


- — — — — 一 3 
S 22 B, Por - Bra Bas m. 

Eqns. (35) and (36) give e and X for any combination of applied 
P* and T* if membrane theory alone adequately defines what is taking 
place. This solution will contain an appreciable error if the bending 
stresses, which the theory assumes to be zero, are of sufficient mag- 
nitude. Any such stresses would have a negligible component in the di- 
rection of the applied tensile load but they may constitute a significant 
portion of the torsional component. If it is assumed that bending 
stresses have a negligible effect on P*, membrane theory alone will 
adequately describe P* as in eqn. (31). If the cross-coupling coeffi- 
cinde, Gii- and CN Au, retain their equality in eqns. (31) and (32), 
then the only significant departure from membrane theory would be 
evident in eX лә" 

Assuming “0 ex 15: and D An are all constants as determined 
by membrane theorv, a different theoretical basis is chosen for 
evaluating 20. 
From eqns. (34) and (36) there is 


ex nu 


дә — . (38) 
б Za ex me = در ی‎ 
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ex 5» can be written as 


2 2 
Xu təx ۔‎ Kia d es 


> وو کی 


a (39) 


X22 "EC X 
22 ox, 


In the theory of torsion for prismatic bars of narrow rectangular 


cross section the angular deformation is related to the applied 


torque m. 


Y | 
= m m “— — —— —-. A 
e Gh'(.bb7a -.21h) (40) 


ҹ 


A recent investigation by Chen chee has modified this relationship by 


a factor, f, so that it applies to rectangular bars with an initial twist. 


ES : an 
T — fGh3i(.667a -.2Ih) l 


where 


Е = | + (1+ 4) ( .زا‎ EE)" (42) 


L mo site YƏ) S,, and J. N, Goodier, Theorv of Elasticitv 
(New York: McGraw-Hill Book Company, 1951) p. 278. 


57 Chu, The Effect of Initial Twist on the Torsional Rigidity 
of Thin Prismatical Bars and Tubular Members, Proceedings of the First 
U. S. National Congress of Applied Mechanics, (New York: The American 
Society of Mechanical Engineers, 1952) p. 265. 
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Non-dimensionalizing eqn. (41) there is 


6 3ه 2 


— —X ^ 
T: С.Һ? (. 667a. - د + اإ(طا2,‎ += Z (144)( Omar = 214 


For the condition of zero applied tensile load in eqn. (36), 


x 
B22 = Tx Es 
Eqn. (43) is therefore an expression for Ê دد‎ based on the theory 
of torsion. Since this relationship is for a rectangular cross section, 
a slight modification of the h term is in order for the present case. 
As a first approximation it is assumed that the average blade thickness 


normal to the z-axis, 77 may be substituted. 


Ы 
Ы 4 - 
0 
haya > h (1-14? :)dg s h(1- ER? C 
avg i $ f ә 3 
0 


Substituting into eqn, (39), @ 22 is found to be 


Chavg (, 670-21 hag) 147% (144 )( Grae =) | uo 
2.a5E - 





X 22 7 


(47) 
This establishes the last unknown influence coefficient in eqns. 
(31) and (32), allowing a complete determination of these equations as 


well as eqns. (35) and (36). These equations now completely describe 
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the deformation of a twisted blade under any combination of applied 
tensile load and torque, within the range of validity of the previously 


listed assumptions. 
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3, Test Equipment and Procedures, 
Purpose 


The purpose of testing was to obtain experimental data for longi- 
tudinal and angular deformations for known values of applied torque and 
tensile load utilizing blades of different initial twist rates. From 
these relationships experimental values were available for: 

a. Young's Modulus, E. 

b. The strain distribution across the chord for comparison with 
the parabolic distribution predicted by eqn. (25). 

C. 6 coefficients in eqns. (35) and (36) for comparison with 


the theoretical values. 
Background 


The basic requirement for testing a twisted blade under the con- 
ditions described in the Analysis is to have a test apparatus which 
allows the application of arbitrarv, but known, external tensile loads 
and moments. Lacking any such apparatus, tests were devised to inves- 
tigate the two extreme situations: 

Test I: An applied axial load with no applied moment. 


Test II: An applied moment with no applied axial load. 
Test Specimens 


Test I can not be accomplished in a normal tensile test since any 
tendency for the blade to untwist in reaction to applied tension would 


be resisted by the gripping head, thus introducing an undesired applied 
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torque of unknown magnitude. Even a 'swivel" head would tend to intro- 
duce such a torque due to the large frictional forces involved. For this 
reason the test specimen was machined to twist in one direction over 
half the length of the specimen and then the twist was reversed over 

the other half. These halves were therefore symmetrical except for the 
reversed twist. Any tendency for one half to untwist under tension was 
exactly equal to that of the other half. The center of the specimen was 
then free to untwist and, since one twisted section had no effect on the 
other, each half was essentially a twisted blade in tension but free of 
applied torque. 

The test specimens were machined with a 1.25 inch diameter cutting 
bit from one and one-half inch diamater 7075-T6 extruded aluminum rod. 
The design was such that a cross section normal to the z-axis was 
essentially a rectangle of width h = .2 inches and length 2a = 1.5 inches. 


Thus ? = R/h = 3.125. 


Three different test specimens were machined, as shown in Fig. 4. 
The machined twist, as characterized by dO/dz, was 22.5, 11.25, and 
zero degrees per inch for specimens #1, #2, and #3 respectively. At the 
center of specimens #1 and #2, prior to reversing the twist, a short 
section of the rod was left intact to allow for the measurement of the 
deformation twist, Y. The twisted sections were approximately nine 
inches long so that strain readings taken at the midsection would be 


free of end effect. 
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Longitudinal Deformation 


Longitudinal deformations were measured by means of strain gages, 
type FAP 12-12 (S13) and type A-7, bonded to the specimen at the center 
of each machined section, Placement was such that strains along the 
tangent to the helix were obtained at $ = 0, 1/2, 3/4, and 1 on each 
side/edge of the blade. This duplicate placement plus the "double- 
ended'' feature of the two twisted specimens allowed for a good 
multiplicitv of data. 

The strain gages were fed through a 20 channel Baldwin-Lima- 
Hamilton Switching and Balancing Unit and into a Budd DATRAN Digital 
Strain Indicator, shown in Fig. 5. This arrangement combined the 
advantages of a direct strain reading in micro inches per inch with a 


rapid reading capabilitv on all gages. 
Angular Deformations 


If a circular disc of radius 11.46 inches were attached as a 
collar to the specimen, a radial angle of five degrees would subtend an 
arc length of one inch on the circumference of the disc. An angle 
indicator, shown in Fig. 6, was devised with this principle in mind. 
Rather than the complete disc, a three inch width of the disc was 
clamped around the specimen and the circumferential arc was calibrated 
in inches. A sighting gage fixed in a suitable stationary position then 
allowed circumferential linear readings for angular deflections. These 
readings of Y , taken in sixteenths of an inch, were converted to A 


by the multiplying factor 
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ela 


L-g n ETY as 


where L is the length of the twisted section, 


Test 1 


Test I was conducted in a 300,000 pound Riehle Testing Machine, 
The specimen was positioned to accept a tensile load and gripped in 
fixed clamping heads at either end. The angle indicator was clamped 
between the two twisted sections for specimens #l and #2, with the 
sighting gage appropriatelv placed to read instantaneous angular 
deformations. The test set-up is shown in Figs. 7 through 9, 

The strain gages and angle indicator were zeroed at a preload of 
1000 pounds, referred to as "gage zero", Readings were taken at 2000 
pound intervals to 10,000 pounds gage, This range of loads allowed 


the maximum strain to remain well below the elastic limit of the material, 
Test 1 


Test II was a torsion test with P = O, For this test a torsional 
apparatus was set up as shown in Figs. 10 through 12, The twisted 
specimens were clamped at the center section and supported at the free 
end by an idler, in which the specimen was free to rotate, Also 
attached at the free end was the angle indicator and a moment arm. The 
moment arm was designed to act as a ten inch lever when loads were 
applied at a right angle to it یا‎ oh the eye bolt. A cable attached to 
the eye bolt fed through a dynamometer and pulley arrangement to a 
take-up ratchet. In this manner the applied moment was simply the 


dynamometer reading multiplied by the ten inch lever arm. 
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As in Test I, an initial load was applied and all gages zeroed. 
Readings of strain, angular deformation, and applied torque were then 


taken in increments of applied torque. 
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4 Results and Discussion. 


The experimental value of Young's Modulus, as determined from 

Fig. 13 and Hooke's Law, was E = 10.1 x 10° psi. This is two percent 
lower than the handbook vaite" for the test material but was considered 
to be the most accurate figure available for that material sample. The 
experimental E and the handbook values for the shear modulus, G and 
Poisson's ration, jii E were used as material constants in all 
calculations. An experimental value for the shear modulus could have 
been determined by a torsional test on the round stock prior to 


machining. From this test the shear modulus could be found from the 


torsional relationship for circular rods, 


T 


Y E 


where J is the polar moment of inertia. The absence of such a test 
was an oversight in the test procedure. 
Table I lists the =@ and C coefficients for each of the 
three specimens, comparing the theoretical values with those determined 
by test. A graphical comparison of theoretical and experimental values 
IES States Department of Defense, Strength of Metal Aircraft 


Elements: Military Handbook 5 (Washington: Government Printing 
Office, 1955) p. 114. 


STD 


EDU O, W., Handbook of Engineering Fundamentals (New York: 
John Wiley & Sons, 1965) p. 12-56, 
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of the various coefficients relative to the initial twist rate is 
plotted in Fig, 1&, It is significant to note the rather good 

degree of correlation attained by a relatively simple theory. The 
cross coupling coefficients, Cx 12 and ex 21 6 12 and 6 217 

vere found to be equal as the analytical solution predicted, Prelim- 
inary investigation made it clear that membrane theory alone produced 
extremely erroneous values for كه‎ 06 The modification introduced in 
the Analysis to obtain به‎ 22 from a theoretical 8 22 based on the 
theory of torsion proved to allow a reasonable prediction of this coef- 
ficient. It should be noted here that an error was introduced in the 
application of Chen Chu's modifying factor for the high initial rate of 
twist since he stated that there was an upper limit in his analysis, 
beyond which higher order effects were significant, L 1) This upper 


limit corresponds to the initial twist of Specimen #2, 


TABLE 1 


TABULATION OF RESULTS 





Specimen #1 Specimen #2 Specimen #3 
K = ,2945 К = 3 К = 0 


















11 , 866 892 , 963 , 947 1.0 . 980 
12 .0749| .0668 .0389 

21 , 0749| 6 : , 0396 

22 , 02041 9 1 ,0115 

m 1.69 | 1.62 | LE 

12 -6.21| -6.40 s 4,16 

21 -6,21| -6.67 a 4,24 

22 71.8 85.5 101.3 
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The three experimental values of 6 11 Were obtained from Fig. 15, 
where e, versus P* is plotted for T* = O (Test I). The other coefficient 
obtainable from Test I was B 217 the ratio of A to P* with T* = 0. 
These values were obtained from Fig. 16. Likewise, Test II provided 
experimental determination for coefficients 6 12 and B 227 


e oa, and (X /T*) y» respectively, as obtained from Figs. 17 and 


Р*= 
18. These B coefficients were converted to experimental ex co- 
efficients by eqn. (37). 

Figs. 19, 20, and 13 are plots of strain data obtained from Test I 
for Specimens #1, #2, and #3 respectively. The strain reading from 
every station, modified by the factor (1 + K” 6 7), vas plotted versus 
g for P = 10,000 pounds, T = O, This strain distribution across the 
chord is seen to plot in a linear fashion as predicted by eqn. (25), 
which is drawn as the solid line. Since the tangential strain at all 
chordwise stations was found to vary linearly with load, the value of 
10,000 pounds was simply a convenient basis and any other load would 
plot qualitatively the same, 

The strain distribution over the chord was essentially parabolic 
for twisted specimens, with the greatest strain realized at the mid- 
point of the chord. The strain relief evident at the outer fibers is 
directly attributable to the amount of angular deformation, or "un- 
winding" of these specimens., An interesting point to note in Fig. 19 
is that while the specimen is subjected to a 10,000 tensile load, the 
outer fibers are actually in compression. This is a result of the extreme 
angular deformation for the case of high initial twist rate, For the un- 
twisted specimen the strain distribution was uniform over the chord, as 


one would expect. 
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Figs. 21 through 23 are similar plots of strain distribution across 
the chord but for the condition T Z 125 ínch pounds, P z O. Again the 
parabolic nature of the strain distribution is evident from the data 
for the twisted specimens, but there is no discernible pattern for the 
untwisted specimen. This last case appeared to be the only test in 
which the measured surface strain was not representative of the actual 
strain at the central fibers at each station. The "back to back" 
arrangement of strain gages on all specimens served as a check for this 
type of discrepancy. Test II of Specimen #3 was the only case in which 
a noticeable difference in strain readings was found from one side of 
the blade to the other. These differences were qualitative as well as 
quantitative, with es actually changing from tension to compression from 
one surface to the other. Even an average of each "back to back" pair 
of strain gages failed to provide a discernible pattern of behavior 
for this specimen. 

The reason for the lack of conformity of this test can most likely 
be traced to the test apparatus. The torsion idler ideally allowed the 
specimen to rotate completely free of friction but, at the same time, 
completely supported in all directions so no bending normal to the 
z-axis would occur. A compromise was adopted for this item in that the 
idler was simply a block support with a slightly oversize hole to allow 
free rotation. The resulting "play" in this hole allowed some degree of 
bending of all specimens as the moment was applied, but only in the case 
of Specimen #1, where the magnitudes of longitudinal strains due to 
torsion were so small, did the longitudinal strains due to bending 


become an important factor, 
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The untwisted specimen, of course, violates the assumption of 
"large" initial twist rates that was made in the analysis of membrane 
stresses. However, at low rates of twist the membrane stresses become 
insignificant in comparison with the shearing stresses in resisting an 
applied moment. The untwisted specimen was therefore included in the 
test as an important limiting case. It also served as a means of 
experimentally determining Young's Modulus. 

The results of these tests have shown the theoretical approach 
herein utilized to be a good first approximation for defining the 
reaction of a twisted blade to a combined stress field. The results are 
not exhaustive, however, since specimens of many different initial twist 
rates and cross sectional aspect ratios would be required to establish a 
definitive pattern of response. An ability to apply completely arbi- 
trary combinations of load and moment in the test procedure would 
provide for a more thorough investigation of experimental-theoretical 
correlation. 

Second order effects which were neglected in this analysis assume 
much more importance at very low initial twist rates. It is assumed 
that membrane theory would prove more reliable for a thinner blade 


cross section; i.e., as the blade actually approached a pure membrane. 
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52 Conclusions. 


The conclusions to be dravn from this paper are, for a tvisted 
blade of narrov rectangular cross section: 

a. Modified membrane theory provides a good first approximation 
to the response of the blade to combined loads of tension and torsion. 

b, The strain distribution across the chord is essentially 


parabolic and is accurately predicted for the tvisted specimens tested. 
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6 Recommendations. 


It is recommended that further work in the area of twisted blades 
under combined loading be directed toward 

a. amore complete theoretical treatment, to include, for 
example, a detailed analysis of the shear stresses in the plane of 
the cross section as well as an investigation of the effect of higher 
order terms that were neglected in the present analysis, 

b. the development of an experimental apparatus that will 
allow a complete freedom of choice of combinations of the applied 
loads, and 

c. a confirmation of the validity of the conclusions made in 
this paper by tests on a number of specimens in which both cross 


sectional aspect ratios and initial twist rates are varied. 
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Fig. 2. The Blade Thickness, t. 
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Fig. 21. Strain Distribution Across the Chord, Specimen #1, 


Test II (T = 120 inch pounds). 
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Fig. 22. Strain Distribution Across the Chord, Specimen #2, 
Test II (T = 129 inch pounds). 
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APPENDIX A 


TABLE IT 


TEST SPECIMEN SPECIFICATIONS 


specimen” Cross Section 


ә 
(in.)|(in.)|(in. 


#۴1 7 ,2945 . 2864 
#2 4196 9:4 . 1963 . 1473 .1473 
13 ,200 12.7 0 0 0 


T = R/h = 3.125 for all specimens 





1 
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TABLE 1 
STRAIN GAGE PLACEMENT 


SPECIMEN #1 


Gage! EE Cage" 
Number Type Factor  |Distance 





— so the "back to back" gages appear as a pair. 


2 "uA" corresponds to type A-7. "B" corresponds to type 


FAP 12-12 (S13). 
3411 gage resistances are 120 t s ohms, 


全 ev gage. 
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TABLE IV 
STRAIN GAGE PLACEMENT 


SPECIMEN #2 





Ha. Gag = Cage” ¢ 
Number Type Factor 6۵ 





O ou LFW N pa 


DD DD GG DD SS DD DG DG DD DS DD D حم‎ 





so the ''back to back'' gages appear as a pair.‏ کو و 


Сид" corresponds to type A-7. "B" corresponds to type 


FAP 12-12 (S13). 


EMT gage resistances are 120 t .5 ohms. 
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TABLE ۷ 
STRAIN GAGE PLACEMENT 


SPECIMEN #3 










as ۱ Gage 
Number | Type 






GES, 
Factor 


t 


Distance 












== 
ON ww‏ من سم سم 


= 


= 
non su 
بن‎ 





Listed so the "back to back" gages appear as a pair. 
- gage resistances are 120 t .5 ohms. 


cy gage. 
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TABLE VI 


TEST I DATA, SPECIMEN #1 


Y (1/16 EEN 










P (pounds) 
P* x 10? 





xX x 10° 


Gage Number س‎ Strain (microinches per inch) 


1 350 715 1080 1460 1860 
2 325 670 1020 1400 1780 
3 1100 2200 3280 4370 5440 
4 1075 2150 3220 4280 5340 
5 280 580 880 1200 1530 
6 350 720 1100 1490 1890 
7 325 650 1000 1350 1720 
8 325 660 1010 1370 1750 
9 1020 2030 3040 4040 5040 

10 1120 2230 3350 4460 5570 
11 300 600 920 1250 1600 
12 310 630 960 1310 1670 
13 860 1730 2600 3490 4380 
14 830 1650 2480 3320 4160 
19 810 1620 2430 3260 4090 
20 750 1520 2290 3080 3880 
25, 820 1640 2490 3350 4210 
26 
31 740 1490 2250 3040 3820 
32 800 1620 2440 3270 4100 
37 -190 -350 -490 -600 -700 
38 -150 -280 -390 -480 -540 
39 -170 -320 -460 -570 -660 
40 -200 -380 -540 -670 -790 
É y 102. 1.08 2.15 3.22 4.29 5.35 


— سس‎ 一 一 — — mg ےہ‎ -: 


Los. eqn. (48). 


Faulty gage. 
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TABLE VII 


TEST I DATA, SPECIMEN #2 





adına 2000 4000 6000 8000 

P* x 10° 0.66 1.32 1.98 2.64 

25.2 19.2 12.8 6.40 !ہو فو می 

X x 102 2.85 5.70 8.54 | 11.21 . 
Strain (microinches per inch) 





1 1610 2160 2720 
2 1680 2250 2820 
3 2440 3245 4060 
4 2450 3265 4080 
5 1570 210 2680 
6 1650 2240 2840 
7 2070 2760 3460 
8 2080 2790 3470 
1980 2660 3345 
2120 2830 3560 
2100 2800 3480 
2120 2830 3550 
2310 3070 3820 
2240 3000 3740 
1680 2250 2830 
1520 2040 2570 
2530 3370 4220 
2350 3140 3930 
1740 2330 2910 
1620 2190 2760 
1235 1670 2120 
1120 . 1540 1970 
1100 1510 1930 
1190 1620 2055 


lee eqn. (48). 
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TABLE VIII 


TEST I DATA, SPECIMEN #3 


— — 


۳۷۷۵۰۵2 ( 2000 4000 6000 8000 İ 10000 
px x 107 N66 1.32 1.98 2.64 3.30 
Y (1/16 inchi 0 0 0 0 0 
xylo” 0 0 0 0 0 


Gage Number Strain (microinches per inch) 





= 


= 
Nn NU RP OA ON © 


1 





2190685 . (489% 


Aene gage. 
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TABLE IX 


TEST II DATA, SPECIMEN #1 








T (inch pounds) 


T* x 109 91.4 
t afis mo 1727 
X, x 107 7.79 


Gage Number Strain (microinches per inch) 





m [o [m [m m fo‏ مہا یں 


Hs eqn. (48). 
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TABLE X 


TEST II DATA, SPECIMEN #2 






T (inch pounds) 24 49 
Ty x 10° 10.56 6 
عو م نب‎ Inch) 
A x 10° 


67.32 |80.08 





ves o [e 0 | m m Poe 





eqn. (48).‏ سز 
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-—Á— HÀ 


107 126 
47.1 | 55.4 
1۳۰ 2 19.1 

6.15 


5 428 


TABLE XI 


TEST II DATA, SPECIMEN 3 


16 33 49 68 86 
7,0 | 44.57) 25560 | 29.9 à 3249 
2.4 5.2 7.5 | 10.4 | 13.4 
.77 | 1,67 0 2842 | 3.35 | 4.31 


Strain (microinches per inch) 


-30 -36 
36 45 
-20 -24 
84 99 
-79 -93 
ell -12 
-38 -44 
93 110 
“147 | -171 
-26 -31 
43 48 





T (inch pounds) 

T* x 102 

Y (1/16 inan 
| X x 102 


Gage Number 







— 


= 
NO Jin meme WH ODM WOW 
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E eqn. (48). 


E gage. 


APPENDIX B 


Calculation of Young's Modulus, E. 
From Fig. 13: 


= 10,000 pounds 





e. = 336 x 1077 inches/inch 
avg 
P P 10000 6 : 
5 “Ae, ~ 2ahe. ^ 2(.75)(.2)(.000336) ^ 10.1 x 10 psi 


Calculation of ek Coefficients. 
The calculation of the first three ex coefficients involves the 
solution of the appropriate integrals in eqns. (29) and (30). These 


can be broken down into basic integrals of the form 


fe, ie +g ә ар 


where n = 0, 2, and 4. These integrals will be referred to as Ij, L, 


and Ij; respectively. 


l 
۱ | | 
ei [eg ав 2 201722 to GR tank 
o 


1 2 4 6 
CSS = 1 - K + K - K + م .امه .امام‎ 


but 


and tan”? K - K - ik + 2 k Кр 


Substituting, 
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B: ۷)۳ e ә ( l-gk + 二 K^] 


I. = Ltrs 


= 


N 


۱ 2 
әр Буд. al 
nee ` Sege дә tan" 
»i 
owe zl Be لے ۔ خم عط + 4 ہے کڈ اب‎ Kt...) 


i 


دا 


ə k^ vv 
4 
CH: | 3 
هيلو سم رد1‎ - — 1 
ə ETS عم‎ > alien) x 
4 
ات - ود‎ )۱ 2+۱۰) + ۷ - tee] 
"ҹа ев бе, 0, 
I. = = > K * 3 K + 
From eqns. (29) and (30), 
cuc le pk I; z (1 22 «۰۰ - 2) ۱۷۶-2 Ли 
دک 2س 3) لط + ۳۶ +2) + - ۱ 2 وروی‎ 
X12 = e a2, 7 log تر و‎ 


۱ 3 
c2 = (3 K - ZK +5 Kre) 0 دم 3 ذم‎ + ,۰( 
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° = y K - # (2+%) K + “x 7 


From eqn, (47), 
ور‎ G h ). "۲7 -.21haq)|[! + (1222) (é. Sas TN. 
22 
2 a ع‎ 


e (47) 


Substituting into the above three expressions for ex 11» ex,» and 


22° there is 


o) = 1 = .148 + .014 = ,866 
c), = 1 - .037 = .963 

p ! 

c5), » .0982 - .0262 4 .0029 - .0749 


046. = 0033. - ۰0491 = )ر1 


) 5,9 X/0 72.002٦ 
27 194000 ر٤۶‎ 


[.cor(.75)-.21 (182) J[ 


[i+ 2 (14.33)L 2004 25) 1. HA 


12 


501 = .01516(.5 - ,0382)(1 + ,988) + .0065 
ox - - 

0204, = 0065. + 0139, = 221 
0022. 4 )222. 4 0412)(1, - 0176(.5. = و( وروی 


0121. = 0022. + 00986. = كه 


ex = .0183(.5 - .042) (1) + 0 = .00838 


و221 


Calculation of B coefficients. 


The expressions for B coefficients in terms of æ coef- 


ficients are to be found in eqns. (33) and (34). Let 


.2 
تر = و رت D - e,‏ 


dv 
Bu = 


2 رک سه 
Pu. =p‏ 
B 22 * D‏ 
D, - (.866)(.0204) - (.0749)” - .0121‏ 


B 1171 2 (.0204)/(.0121) = 1.69 
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B (1 = (-.0749)/(.0121) = -6.21 


= (.866)/(.0121) =. 71.8 


“o 
N 
N 

— 
سم 
۱ 


О, — (.963)(.0121) = (.046) = .00954 


= (.0121)/(.00954) = 1.27 


Wie 
نس‎ 
F= 
< 
N 
۱ 


2۔ 


Br» = (-.046)/(.00954) 


= (.963)/(.00954) = 101 


US 
n 
N 

١ص‎ 
NO 
۱ 


7 = (1) (.00838) = .00838 


(.00838)/(.00838) = 1.0 


"DO 
= 
= 
“— 
дә 
II 


(.963)/(.00838) = 119,3 


B 2273 


Calculation of theoretical strain distribution. 


The theoretical strain distribution across the chord is given 


by eqn. (25). 
e. (1 + Kg") = e 4 XK“ (25) 


14 


a) Test I: P = 10,000 pounds, T = O 


From eqns. (35) and (36) 
“ə” Pu Dy 
X= Ba P* 

Specimen #1: 


1.69(3300 x 107°) 455593 m agi? 


M 
ә 


-6.21(3300 x 1079) = 720500 x 107° 


۷ 


X 


l 


e (1 + Kg") = 5577 x 1079 - 20500 x 1075(.2945) g^ 


6 2 


22 E -6 
e.(1 + Kg") = 5577 x 10” - 6035 x 10 £ 


Specimen #2: 


1.27(3300 x 1079) 2 4191 ser” 


0 
ib 


-4.82(3300 x 109) = -15,906 х 1072 


X 


e (l+ KE) s 4191 x 107” - 15,906 x 1077(.1473) g” 


2.2 ар 2 
e.(1 + K€") = 4191 x 10 - 2343 x 10 f 
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Specimen #3: 


1.0(3300 x 107°) = 3300 x 107° 


0 
۱ 


DR -6 
e, (1 4 K'g^) - 3300 x 10 


Test II: T % 125 inch pounds, P « O 


From eqns. (35) and (36) 


io B2 Es 


Specimen #1: 


6 


O 
۱ 


-6,212.2 x 10 ( - - 2 Ло 


6 


il 


X = 71.8(52.2 x из On 


e.(1 + KE") = -324 x 107° + 3748 x 1077(.2945) g^ 


6 


202 - -6 م‎ 2 
e, (1 + K £ ) = -324 x 10 + 1103 x 0 f 
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b) 


Specimen #2: 


-4.82(56.76 x 107°) ددویے‎ یو٤‎ 


(D 
li 


104 رر یرت 18 66.76( TON‏ 


x 


6 


e (Ë + KE) 0 0 08 107”(.1473) £” 


-273 x 107” 


22 -6 2 
e,(l + KF ) + 844 x 10 م‎ 


Specimen #3: 


6 


X = 119.3(55.4 x 107) = 6610 x 10. 


—Ó‏ مسر 
)F = О‏ 10 * 066610 ع e (l + KE)‏ 
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